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AiisrmcT 


A tlieoretical study of tiie linear stability 
(both temporal and sxoatial) of flow in a rigid, concen- 
tric annulus to axis ymi^ie trie and non-axisymmetric dis- 
turbances is presented here. Tlie continuity and Navier - 
Stokes equations are reduced to a set of coupled, linear, 
ordinary differential equations goveriring the propagation 
of the disturbance. The flow is said to be temporally 
stable if the disturbance, applied at an initial instant 
cveryidicre in the fluid and periodic in the downstream 
direction, decays with time while it is considered to be 
spatially stable if the disturbance, periodic in time 
and imposed at a specified location, decays in the down- 
stream direction. The coupled differential equations are 
solved numerically for the complex eigenvalues . The in- 
tegration starts at the outer wall by means of a Frobenius 
scries and is continued to the inner wall by the fourth 
order Eunge-Lutta method. An eigenvalue search technique 
is employed to ascertain the number of oigenvaEues within 
a closed region of the complex w or k-plane. Actual 
isolation of the eigenvalues is achieved by means of an 
iterative scheme. Eesults are presented for the linear 
spatial stability of the annular flow to axisyQimotric 
disturbances. Throe different annuli with diameter ra- 
tio y = 0.2, 0.5 and 0.8 are considered. It is found 
that the flow is spatially stable to infinitesimal axi-r 

i' 

symmetric disturbances up to a modified Reynolds number 
of 10,000. 



CHAPTER 1 


INTRODUCTION 

Numerous efforts have "been made to clarifj^ and 
explain theoretically the reinarkahle process of transition 
from laminar to turbulent flow. The recent theoretical 
investigations are based on the assumption that lamina.r 
flows are affected by certain small disturbances; for 
example, in the case of a boundary layer on a solid-body 
placed in a stream, these disturbances may be due to wall 
roughness and/or due to irregularities in the external 
flow stream, whereas in the co.se of pipe flow, they may 
originate in the inlet section. The sto.bility theory en- 
deavours to follow the behaviour of suc.h disturbances, 
when they are superimposed on the main flow, the decisive 
question to answer in this respect is whether the distur- 
bances amplify or decay. If the disturbance ultima tel 3 ?' 
decays to zero the flow is said to be stable, but if the 
disturbance amplifies the flow is said to be unstable. It 
does not follow, however, that instability leads to turbu- 
lent motion. 

The behaviour of small disturbances to a flow 
can be analyzed in two different \TB.ys , In the first of 
these, a disturbance is assumed to be applied at an ini- 
tial instant everywhere in the fluid and is, in particular. 
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periodic in the downstream direction. Tlie flow is then 
said to he 'temporally’ stable if the disturbance decays 
with time and unsta,ble if it grows. In the second appro- 
ach, it is assumed that a disturbance which is periodic 
in time is imposed at a specified location in the fluid 
and propagates doxms tr earn . If the disturbance decays in 
the down-stream direction, the flox^ is considered to be 
'spatially' stable. Though the latter point of view is 
physically more realistic, alsiost a3.1 effort has been devo- 
ted to the former apj)roach since the very inception of the 
stability theory. A possible reason for this bias is given 

at the end of this chapter. lAlhile the stability of flow in 

❖ 

a circu].ar pipe has been analysed extensively [ 1 - 5 ] 5 
much atcontion has been paid to the flow in an. annulus. 

Only the temporal stability of flow in a concentric annulus 
to infinitesimal axisynmetric disturbances was considered 
by Mott and Joseph [o ] . 

The purpose of this study is to analyze theore- 
tically the linear stability of viscous, incompressible 
flow between coaxial, right-circular cylinders, Both 
axisymmetric and non-axisytmiietric disturbances are con- 
sidered. 

*Numbers in square parentheses designate references listed 
at the end of the dissertation. 



The corresponding mathematical proTolem can be 
expressed as follows ; 

For time t < 0 we suppose that the flow is 
undisturbed; for t > 0 one of the following two situations 
exists : 


(i) For temporal stability'-, a disturbance is assumed 
to be applied everjntiere in the fluid at t = 0. This dis- 
turbance is periodic in the dovmstream direction z and 
is an a.rbitrarj'" function of radial and angular coordinates 
r and 0 . 

(ii) For spatial stability, a periodic disturbance is 
imposed at z = 0 and at t = 0. This disturbance has a 
frequency cu and is an arbitrary function of r and 0. 


The 0- dependence of the disturbance velocity 

3|c 

vector V can be expressed as a Fourier series in 0 of tho 
form 


V * (r,©, z, t) = Re 


V (r, 0,z, t), j = Re 


oo 

z: 

‘'n=0 


Vn (r,z,t) 


^in© 

® J’ 
( 1 . 1 ) 


where n is the azimuthal mode number, Re denotes the 
real part of a complex function, and i = v/-T . For n = 0, 
we have an axisymmetric disturbance ; for n = 1 , the per- 
turbations are constant along a simple helix; foi* n = 2, 
the perturbations are constant along tv/o inteivined 
helices, and so on. 
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For tlie (r,z,t) dependence, there are two 
possible \m.js of repres'enting (r,z,t) depending upon 
the tj^e of stahility analysis. 

(i) For temporal stahility analysis, each (r,2,t) 
can he expressed as a Fourier series of the form 


oo 


(r,z,t) - ^ 2;^^ (r,2) e-^Hi^ 


m =1 


^nm 


(1.2) 


where is the complex frequency to be found. Finally, 


the z- dependence of v (t,z) C8.n be expressed as 

liil 

V (r,z) = V (r) e^^^ , 

■"nm ‘ — nm ^ ^ j 


o-ron 


( 1 . 3 ) 


where k (an imaginary quantity) is the vrave number of the 
disturbance. Thus, the disturbance velocity vector is 
given by the double Fourier series as 


r c>o CO 


v*(r,e,z,t) = Re Znm ^ 

( 1 . 4 ) 

(ii) For spatial stability analysis, each v (r,2,t) 
can be expressed as a Fourier scries of the form 


00 

Zn <^^>2,t) = X 


m =1 


Mill 


k z 

,, (r,t) o , 


( 1 . 5 ) 


where is the complex i.\ra.vc number to be found. Finally, 
the t-dencndencc of v^^ (r,t) can be expressed as 

^liin ' 


^ ^ j N = / N -i w t 


-nm ’ ~nm 
where tJ is the real frequency of the disturbance. Thus, 


( 1 . 6 ) 
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the disturbance velocity vector is given by the double 
Fourier series as 

r 5*^ oo ^ u 

v^Cr.OjZjt) = Re XI J ^ • 

(1.7) 

Note that V (r,9,2,t), v^ (r,z,t), v^ (r,z), v^ (r), 
v,, . . , (r,t) and v^._ (r) arc all complex vectors. 

The flou is said to be temporally stable if the 
imaginary part of the complex frequency while it is 

said to bo spatially stable if the real part of the complex 
■vra.ve number k < 0. Another Important difference bctv/cen 
the two cases is that, because of the continuity equation 
(soe equation (2.20)), the vector eigenfunction^ (r) is 
equivalent to two scalar eigenfunctions for the temporal 
stability rather than three as in the case of spatial sta- 
bility. The problem of spatial or convective stability is, 
therefore, much more complicated than that of temporal 
stability, which is quite likely the reason why it has been 
neglected in the past in favoiir of temporal studies. 
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THEOEETICAL ANALYSIS 

2 . 1 ]yfathematlc3.1 Formulation 

It has already heen pointed out that the stahllity 
of a flow is determined hy assuming small disturbances 
superimposed on the basic flow. Therefore , denoting the 
basic flov'j" velocity vector by ^ and the disturbance velo- 
city vector by v* , we can write the total velocity 
vector \ as 

x=i* + s*, (2.1) 

and similarly the total pressure P as 

p = p* + p'"^ (2.2) 

Since we are considering the sta.bility of laminar 
flow between tvro coaxial, right-circular cylinders, the 
basic flow has only an axial component of velocity and this 
component is only a function of the inner and outer radii 
of the annulus, i.e. J, = (r) Other principal 

assumptions for the analysis to follow are 

(i) Incompressible, New'tonian fluid with constant 
coefficient of vi.scosit 3 ^, 

(ii) No body forces, 

(iii) Steady basic flo^ir, 

(iv) Infinitesiml disturbances. 






(v) Rigid, Impermeable cylinders, 
(vi) No slip at the cylinder imlls 


For such a fluid, the equations of motion in 


cylindrical 'coordinates (figure 2.1) are 


Continuity t 


Sr r r SO Jz " " 

and momentum (Navier-Stokes equations) 

( ^ + V + — 2 -:rS - + 


( 2 . 3 ) 


+ V • 

^ t r 


"^9 

r 

^9 ~ 

Y ah ^ 1 

ar 


'i>\ 

Y 

r 

+ i- 

f h 2 ah 

'^02 r2 

2 

a 

^r.2 

0 2 

UiF - 


r2 

^ 2 
r 






(2.5+) 

-4^ + V 
it r 

^r 

’^o 

r 

^V9 

ae 

, h h , , ah 

r 2 ^ z 

- ) + 

^ Pr 


R “IF 

(2.?) 


( 2lz ^ ^ Jq ^ ^ ) -f 1 ^ 

^ •‘&r r B0 ^ f ^2 


^ R , 1 , 1 


1 a R , a ^z 


^ + i; — ^ ) 

^92 ^^2 ^ 


( 2 . 6 ) 


where P is the density and U is the kinematic viscosity 
of the fluid. 



'9 


The steady hasic floi; itself satisfies the e(3_uations 


dr 


^ — * 
^P ^ 

Q 


0 


and 


1 ^ . 1 

? dz ^ ^2 r 


P* (r, 0, z) = P*' (z)j 




f (2.7) 


Br 


so that the basic floi; has the well-kno'v/n form 


V 

z 

2 

[(r§ - rt in ^ + 

{ rf 

- r2 ) in ^1] 

ai 


tv 

1 

o o 

(r 2 + rp In ^ - ( 

1 

C\i CM 

U 


h 

1* = 

r 


= 0 , 



k2.8) 


where ‘T „ is the average, over the flow cross-section, 

O- V '* ' 

— 

of the axial velocity , and r^ and are the inner 
and outer radii of the anniuLus respectively. 


If we now substitute equations (2.1) and (2.2) 
into equations (2.3) through (2.6), nBke use of equation 
(2.7), and neglect the product terns of disturbance velo- 
cities with themselves and v/ith their spatial derivatives, 
we obtain the following set of linearized hydrodynamic 
equations for disturbance velocity and pressure 
* 




V 


Sr 

* 


+ V 


JC + 1 
'r r 

-9^ 


“se 


+ 




z _ 


♦ 

Bz 


1 SP 


dz 


(2.9) 


St 

2 * 

- 

Sr^ r 


f hr 






V. 




t r 


+ 


r 


^9' 


w 2 

2 

• + ::: ), 


,2 S 0 


Sz 


2 


( 2 . 10 ) 


0 , 
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❖ 

^^0 


at 


^ “fe 


2 ^ >i« 

= i, ( 

^ ^j,2 * r -Sr 


^ 2 

lo + !_ 1_I' 

r^ r^ 0^ 


5[« 

0 _L 2 


+ 


+ 


2 * 
Tr 

0 


a ■^' 


■ az2 


), 


( 2 . 11 ) 


and 




z . 

'at "^r 




1* 

z 


^ + Y 

ar z 


❖ 

: 5 V 

o z 


1 ^'0 

d z P S z 


+ 


2 * >!' 
, 2i -v. . 

y ( + 1 

^2 r a r 


+ — 


2 =!= 
1- l_Iz 


2 

^ ^z 

^z2 


). 


( 2 . 12 ) 


Equations (2.9) through (2.12) can he made non- 
dimensional hy choosing as the respective units of length, 
velocity and density the parameters r^, V „ and f . Thus let 

d . ' 3. V * 


V. 






V, 

"" V 


av 

0 


av 




L 

V 


av 


^2 ’ 


^ r2 ’ 


t V 


av 




Y_ 


Y = 
z 


P 


Y 5 
av 


P 


f y; 


2 ’ 
av 


\v ^2 


J 


h ( 2 . 13 ) 


The modified Heynolds numher R is related to 


the actual Reynolds mamher R„ (= 2 (r^ - r.)/l>) 

ci a V £i I 

hy Rg^ = 2R (1 -Y)^ -where jf = r^/r^. The dimensionless ba- 
sic flow velocity Y (from equation 2.8) is given hy 

4/ 



11 


= ^ ( 1 - In r' , 

where <<= 2 (1 - / [(1 + In y +. (1 - Y^) ] 

and /3 = 2 In Y / [(1 + JT^) In T + (1 - V^) J . 

Making suhstitutions from equations (2.13) and 
(2.14) into equations (2.9) throu^ (2.12), dividing equa- 
tion (2,9) ky ( V ^yr„ ) and equations (2.10) through 

S, V iZ. 

p 

(2.12) hy ( V^y^r^ ), and dropping the primes from r, z 
and t for convenience, we get the following set of 
linear, partial differential equations for the fluctua- 
tion properties in non-dimensional form 




, Yr ^ 1 iP 

at z az r BG 

1 r * 1 fe 4 . 1- ^Lle 4 . 2_ S 4 . Ize-, 

= H ' 3^2 r 3r - ^2 ^2 3 g2 ^2 39 3^2 

( 2 . 17 ) 

and 


3^2 Jh 4 . rr JZs 4 . ^ 

3 t n Sr h 32 . 32 

= 1 rtlz ^ I ^ + 1. jtZz. 

^ 3r^ r 3 r ^.2 3 q2 3^2 


( 2 . 18 ) 
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Since the equa.tions (2.15) through (2.18) are 
linear, the various Fourier comioonents forming the 
arbitrary disturbance according to equation (1..^) or (1.7) 
8.re independent of one another. Vie can, therefore, con- 
sider only one Fourier component at a time . Thus we hope 
to find solutions of the type 

0 (r, 9, z, t) = Rg 0 (r, 9, z, t) ] 

= Re 1^0 (r) exp (kz + inG « iut)J. (2.19) 

Equation (2.19) is equivalent to four equations by repla- 
cing 0 by v^, Vq, v^ and p successivel37-. All variables in 
the equation are dimensionless, and the superscript ('^) 
indicates a preliminary complex solution which will lead 
to the real solution 0. The quantity 0 (r) is a complex 
function of r alone. Substituting the disturbance 
pressure and velocity components in the form of equation 
(2.19) into equations (2.15) through (2.18) and elimina- 
ting the exponential factor, we get, after some rearrange- 
ment, the following set of linear, homogeneous, ordinary 
differential equations 

(D + ^) V3, f (iVg) + kv^ = 0 , (2.20) 

[d^ + I D ^ ^ aLl-l - + R (kV^ - it*> )}] 


( 2 . 21 ) 
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[ d2 ^ 1 D - { 


n + 1 


- + 


R (kV^ 



(iTg) 


2n - , n xj- 

^ Xp + - Rp 


0 


and 

[ ^ D - k^ + R (kV^ - i ck, )J 

- Rv^ DV^ - Rkp = C , 
■where D is the operator d/dr. 


V 

z 


( 2 . 22 ) 


(2.23) 


For any type of stability analysis the solution 
of the abo've set yields ' 'the ampli'tude functions (r), 

Vq (r)j (r) and p (r) given the ■value of parameters n 
and R. V/hile considering the temporal stability, an imagi- 
nary ■value of k is assumed to find the complex uj 'while in 
the case of spatial stability, a real ■value of w is assumed 
to find the comp3.ex k. This is done by applying the boun- 
dary conditions at the inner and outer walls of the annulus 
to the solutions of the differential equations (2.20) through 
(2.23). It may be pointed out that the dimensional co'unter- 
parts of UJ and k are 


00 * ::: - .. 3 -^..^ 


and 


* 

k = 



( 2 . 2 >+) 


2.2 Boundary Conditions 


The physical restrictions at the rigid, imper- 
meable inner (r = V ) and outer (r = 1) walls are 



and 


( Y ) = Vq ( Y ) = ( V ) = O 5 p ( Y ) finite, 

( 2 . 25 ) 

( 1 ) = Vq ( 1 ) = ( 1 ) = 0 , p ( 1 ) finite. ( 2 . 26 ) 

For given values of real u) (for spatial stability) 
or imagina-rjr k (for temporal stability) and parameters R 
and n, the solution of equations (2.20) through (2.23) 
together v/ith the boundary conditions in equations (2.25) 
and ( 2 . 26 ) leads to an eigenvalue problem. 
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•SOLUTION .OF STAB^ELITY EQUATIONS 

Equations (2.20) through (2.23) together id. th 
the boundary conditions in equations (2.25) and (2.26) 
completely govern the propagation of an arbitrary, infi- 
nitesimal disturbance through an incompressible, viscous, 
coaxial annular flow. Their solution will determine 
whether such a disturbance grows or decays with the doim- 
stream distance z (for spatial analysis) or with time t 
(for temporal analysis). In terms of k or co , the flov; 
would be spatially unstable when k^ '> 0 and temporally 
unstable when > 0, where k^ and are the real 

and imaginary parts of k and u» respectively. 

The set of differential equations (2.20) through 
(2.23) are too complex to be amenable to an analytical 
solution and, therefore, numerical technique has been used 
for the solution. Direct integration of these differential 
equations by means of, say the Runge-Kutta method, is not 
possible since' all the starting values for the eigenfunc- ' 
tions and their derivatives are not known. In order to 
get these starting values, the eigenfunctions are expanded 
as a power series in (1 - r) near the outer wall for small 
(1 - r). The preference of (1 - r) over r for the 
series solution is based on the form of the basic fl'ow 
velocity which being in terms of In r, can be expan- 

ded easily in terms of (1 - p), as 
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(r) = (1-r) + (1-r)^ + ... + + ..., (3.1) 

where = «(+ 2.^ , = (</2) -jS, 

Cjjj = ®(/m 5 m = 3, . 

p 

The expressions 1/r and 1/r in equations (2.20) through 
(2.23) are written as 


and 


1 

r 


1 

,2 


3 ~ r~( l ~r) = ■’ + (1-r)^ + 


+ (1-r) 


m 


(3.2) 


C1 - (1-r)] 


^ = 1 + 2(1-r) + 3(1-r) + 


+ m(1-r)“'"'’+ ... 

respectively/-. Starting from the outer mil [(1-r) = Oj , 
the series solution is carried up to a small but finite -value 
of (1-r) (for example (1-r) - 0.05 i.e, r = 0.95) due to 
practical difficulties of summing a very large number of 
terms. The solution is then continued by- the fourth order 
R-unge-Kut-ta method to the inner wall (r =Y), The whole 
procedure is iterated until the boundary conditions at the 
inner mil are satisfied, thus resulting in a value for the 
eigenvalue k (for spatial stability) or eigenvalue CJ (for 
temporal s tability) , 

3.1 Series Expansion 


Let the power series expansion of v^, 
p be of the form 

Vr = ( 1 -r)^ ( 1 -r) + 7^ (l-r)^ + • • ( 1 -r)“^*’V.. , 

(3»3) 

iv^- ( 1 -r)^ [W^ + ( 1 -r) + W 3 (l^r)^ + ...+\ ( 1 -r)“^‘'‘*+ ... ], 

(3.^) 



% 


17 


V 


z 


(1-r)° 



+ (1-r) + (1-r)^ + 


(3.5) 


and 

p = (1-r)^ fp^ + (1-r) + P^ (1-r)^ + ***■*■ ^ni d-^)^”'^ + • 

(3.6) 


where all V, W, U and P's are complex constants. 

Substituting equations (3.1) through (3.6) into 
differential equations (2,20) through (2.23), we get 
from equation (2.20) 


- (l-r)®'“‘^ aV^ 

- (a+2) ]■*• 




(1-r)^ [v^ - 


(a+1) 7p 

+ d-r)®-”^^ ■ 


7 m-1 

u 



(a+m-1) 7 

J m 


+ V. 


+ * • • 


+ Cl^r)^ nW^ + d-r)''^'^'^ n (W^+W^) + ... 

+ d-r)^'^“^"^ [n{;^ + ••• + (1-:c*)^ kU^ 

^ *1 

+ d-r)^'^'' kU ^ + ... + d - r )'''^^"2 kUj^^^ + ... - 0, ( 3 - 7 ) 

from equation (2.21) 


a(a+i) - aV^ j 
+ d-r)®" [^(a+1)(a+2) - aV^ - (a+l) ” (n^+1) 7^ 

■k (k^+iujR) 7^ j+- d-r)^'^'’ [(a+2)(a+3) \ - a7^ 

-(a+1) 7^ * (a+2) 7^ - (n^+l) (7^+27^) .. RkC^7^ 

+ (k^+iwR) ^2 J"*” +■ d-r)®''‘‘“^~3 ^(a+m-2) (a+m-1 ) 7^ 
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f m-1 ^ o r ^ 

(a+j-1) (n^+1){^ (m~j-1) Y. | 

J J 1 


Kk 


/* ”^-3 ^ 5 ' 

{ ^ Cj Vj-2}-^ V: 


- (1-r)^ 2 nW^ - d-r)’^'^’' 2n (2 M^+W^) 


- (1-r)^'"'”-'' 2n 


m 1 

ZZ (m-3+1 ): . 

0=1 . 

+ (1-r)^-'’ RdP^ + (1-r)^ R (d+1) + ... 

+ R (d+m-1) P^ + . . . = 0 , (3.8) 

from equation (2.22) 

(l-r)^"^ [(b-1) TdW^ j+ (l-r)^""' [l^-CT^+l > -^2 " ^^1 ] 

+ (1-r)^ [(b+1)(b+2) - bW^ - (b+1.) - (n^+1) 

+ (k^+ia>R) wj+ (1-r)^'' [(b+2)(b+3) W^~bVJ^ - (b+1) 

- (b+2) - (n^+1)(W2+2W^) - RkC^W^ + (k^ + iwR) ^^ 2 ]+ . 

+ (1-r)'^'^^''^ r (b+m--2) (b+m-1 )VJ' - T ZZI (b+j-l ) VJ.J 
^ r m-2 . t m»-3 -t 

- “j}- V 3 - 2 } 

+ (k^+luR) \_ 2 l+ ••• - 2 nV^ - (1-r)®-*'' 

m 

2n (2V^ + V^) ~ ... - (l-r)®'"^”^"'^ 2n|^ (m-j+l ) 


+ (1-r)^ nRP., + (1-r)'^'^^ nR (P^ ^2^ **• 


•+ nR ( P.) + ... = 0 , 

-1— • “f ♦i 


( 3 . 9 ) 
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and from equation (2«23) 

(1-.r)C-2 ^(c-1) cU^]+ d-r)^""^ j^c(c+1) - cU^ j 

+ (1-r)° j^(c+1)(c+2) - cU^ " (c+1) - n%^ 

+ (k^+iu)E) U-i ] + (1-r)^'*'"’ [(c+2)(c+3) Uj^-cU^ -( c+1 ) 

- (c+2) - n^ (U^ + 2U^) - EkC^U^ + (k^+iu)R) 112]+ 

r f .^"i 1 

[ (c+m-2)(c+m-1) (c+j-D } 

3 1 

f •» , m-3 

°3 V3-2} 

+ (k^+lu>E) %_2]+ ••• + (1-r)^ RC^V^ + (1-r)^'^'' . 

B (202^ + C^Vg) + ... + 3C3 


b_2 


- (1-r)^ EkP.j - (1-r)^'^'' EkP2 - (1-r)^^^ EkP^ - ... 


d+2 


(1-r)^'^”^"'’^ Ekp^ - ... = 0. 


( 3 . 10 ) 


In equations (3,7) tliTOugh (3.10) any term with either a 
zero or negative subscript' is set to zero. 

To find relations for a, h, c and d, the coeffi- 
cient of r® is set to zero for all s in equations (3.7) 
through (3.10). This gives 

A From equation (3.7) 

either one or more of the three factors (a-1), h and 
c is the smallest, 

(i) If (a-1) is the smallest, 

3.7^ = 0 i.e. either a = 0 or V.j = 


0 
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(ii) If b is the smallest, 

nW^ =0 i.e. = 0 for n ^ 0. 

(iii) If c is the smallest, 

kU^ = 0 i.e. = 0 ( *.* k 0) . 

(iv) If a.ll or any two of these three factors are equal 
but 7^5 and are independent, we get preciselj^ 
the same conclusions as above; in case the corres- 
ponding and are dependent, no useful 

relation can be found for a, b and c. 

B From equation (3 ‘8) 

either one or more of the three factors (a -2), b and 


(d-1) 

is the smallest. 


(i) 

If (a-2) is the smallest. 



(a-1 ) aV^ =0 i.e. either a = 0 or 1 or = 0. 

(ii) 

If b is the smallest, 



2 nV7^ = 0 i.e. = 0 for n 0 . 


(iii) 

If (d-1) is the smallest. 



FidP^ = 0 i.-e. either d = 0 or = 

0 (vR ^ 0) 

(iv) 

Remarks similar to those in A (iv) 

hold if two 


or all of these three factors are equal. 

C From equation ( 3 * 9 ) 

either one or more of the three factors (b-2), a and 
d is the smallest. 
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(i) If (1-2) is the smallest, 

Cb-1 ) = 0 i.e. either h = 0 or 1 pr W^. = o, 

(ii) If a is the smallest, 

2 nV^ = 0 i.e^ = 0 for n ^ 0. 

(iii) If d is the smallest, 

nRP^ =0 i.e. = 0 for n 0. 

(iv) Remarks siinilar to those in A (iv) hold if two 
or all of these three factors are equal. 

D Prom e qua t ion (3*10) 

either one or more of the three factors (c-2), a and 
d is the smallest. 

(i) If (c-2) is the smallest, 

(c-1) cU^ = 0 i.e. either c = 0 or 1 or = 0. 

(ii) If a is the smallest 

RC^V^ =0 i.e. = 0 ( */ ^ O) 

(iii) If d is the smallest 

RkP^ =0 i.e. P^ = 0. 

(iv) Remarks similar to those in A (iv) hold if two or 
all of these three factors are equal. 

Looking over these conclusions, we take 
a = h= c = 1, d=0 and = 0 (3.11) 

so that the boundary conditions at r = 1 are satisfied. 
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The recurrence relations for the constants in 
the series expansions can now he cletcniiined. by substitu- 
ting equation ( 3 . 11 ) into equations ( 3 - 7 ) through ( 3 . 10 ). 

We arrange the terms in ascending powers of ( 1 -r), and 

til 

equate the coefficient for the (m- 1 ) tern in each equa- 
tion to zero to get 

from equation ( 3 . 7 ) 

rr, 1 f p -1 1 r 1 

^ V.^ . n{^ = 0, (3.12) 

J £1 J * 


from equation ( 3 . 8 ) 


( 1 -r) 


m-2 . 


(m-1) mV 


Bli 


{ 


m 


m-4 


m -1 , r ni -2 -i 

3 V 3 }- (-3-1) q} 


J-1 
m -2 

>L 


C. V 
D m 


ni -2 

j=2 


- 3 - 2 } 


.+ (k + itaR) V 


ra -2 


r m-a 

-2n I ^ (m-j-1) W.J-h R(m-1) =0, (3.13) 

3 1 

from' e qua t ion (3.9) 

(^_^)m-2 . ^ (m-j-1) W^.} 


f ^"3 7 P 

^3 Vj-2}-^ (k2+i<*>R) 


m-3 
3 

m -2 
3 


- 2n/lZ (m-o-1) V U nR ( ZZ P.) = 0 , 

L ^^2 V 0-1 ^ 

(3.1^) 
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and from equation (3.10) 


(1-r) 


m-2 


(ni- 


-1) ta-j-DUj} 


r/ 1-2 

I 


-nk{g c. 


u. 


j=i 

ra-3 

I 

5 


ra-2 


^3 Vo-i}-®^ Vi = o> (3.15) 

where again any terra with either a zero or negative 
subscript is set to zero by convention. 

The equations (3.12) through (3*15) are identi- 
cally true for m =1. Since = 0, only and are 

Independent. For determination of all the higher subs- 
cripted Vj U, P and W’s, we note tliat equation (3*12) gives 
°' 3 '(^ation (3.1^) gives W^, equation ( 3 * 15 ) gives U^, 
and equation ( 3 * 13 ) after substituting for from equation 

(3*12) gives P^. Thus replacing m by (m+1 ) in equations 
(3 *12) through ( 3 * 15 )} we get 

.. m 


V 


Y\ ill 

m+1 = TiU7 ( Tmar ( ^ tjItj 


-t m ^ m 

^m+1 ~ rm+lT ^ " ffiTm+i) ^ 


( 3 . 16 ) 




m-1 

j 

m-1 

1 

d 


¥ 


m-1 


iSu7 hh i 5? °3 Vj-1 } 


(3.17) 



u 


m+1 


-} _ j,2 ^ 1 

mIm+TJ ^ m(m+l") t 4^ (j-^J-o) J 


..(k^+i<AjR ) „ R .p .r } 

m(m+1 ) rn-l in(iB+1 ) t / 

J * 


D= 
111-2 


kR 

m(m+1 )' 


m-2 


{ J~ C. U . . 1 
I 0 m-j-l J 


3=1 


+ kR p 
m( m+1 ) 


and 


n 


m 


,2' m-1 




( 3 . 18 ) 


If f “"3 , 


with = 0. 


( 3 . 19 ) 

( 3 . 20 ) 


In that order, equations ( 3 . 16 ) through ( 3 . 19 ) give 
the coefficients V^, W^, and for m = 2, 3, ... in 

terms of and . With all coefficients in the series 

solutions (equations (3.3) through (3 .6)) expressed in terms 
of U^, P^ and W^ , we can express any eigenfunction as a sum 
of three terms 5 for example the radial velocity eigenfunction 
v^ (r) may he written as 

v^ (r) = v^^ (r) + Vy 2 + '^r 3 (^) • (3.21) 


A similar expansion holds for the other eigenfunctions 



3.2 Ste-p-by-SteiD Inte,g:ratlon 
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starting from the outer vjall of the annulus 
(i.e, p = 1) the series expansions developed above are 
carried only up to a small value of (1 _ p) . The solution 
is further continued by a step-by-step integration tech- 
nique to the inner wall ol' the annulus . The equations 

(2.20) through (2.23) are not directly useful for the step- 
by-step integra-tion technique due to the fact that equation 

(2.21 ) contains both the highest degree derivatives in the 
eigenfunction v^ (r) and p (r). To circumvent this diffi- 
culty, we differentiate the continuity equation (equation 
(2.20)) with respect to r and solve for v . The equa- 
tion (2.21) can be solved for the unknown Dp using the 
relation just obtained for D^ v^, and equations (2.22) and 

(2.23) give D (i Vq) and D v^ respectively to yield the 
following set 


2 — 

D V 


( ^ - D) { + n (ivg)} 


1 ' 1 

r 


k Dv - 

z ’ 


(3.22) 


D^ (iVg) = 5 (“ - Rp) + I 


n^^-i-l D , 2 

O ^ - K 


+ R (kV^ - iw )J (iVg), (3.23) 


V 


R (kp + v„ DV ) + / 

T Z L 


£ 2 

2 r “ 


+ R (kf^ - iu> )j. v^ , (3.2lf) 



and 
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Dp = i i D) (IVq) - { 




+ R (k-?2 - iu) )J v.,]o (3*25) 

The integration of these differential equations 
is carried out in terms of the three independent solutions 
of the eigenfunctions, that is, in terms of (r), 
and eigenfunction v^(r), and similarly for 

VgCr), andp(r). 

The fourth order Eunge-Kutta method, used for it, 
is sunmiarized in A.ppendix I. 

3 -3 Satisfying Boundary Conditions at the Inner Wall 


Having integrated the stahility equations over 
the v;hole region V ^ r^ ^ 1 such that boundary conditions 
at the outer vra.ll of the annulus are satisfied (by means of 
the series expansion), it is non required to satisfy the 
boundary conditions at the inner -wall. For n = 0 the 
eigenfunction VgCr) in the equations (2.20) through (2,23) 
decouples itself completely from others. V/ithin the context 
of stability theory, there is, therefore, no production of 
yQ(r) and ue can arbitrarily assume tliat it vanishes under 
the effect of viscosity [7 , page 229}* Thus the set of 
four equations (equations (2.20) through (2,23)) reduces 
to one of three equations only. Thus, the case of axi- 
Asyimietric disturbances differs from that of the non-axi- 
%yrai3<? ■^ric disturbances in the following manner; 
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(i) There are three eigenfunctions andj tliereforej 
three differential equations for the axisyimiie- 
tric disturbance as coirpared to four for the 
non-axisynmie tri c dis turbance . 

(ii) There are two independent solutions for each 
eigenfunction and tuo boundary conditions at 
the inner v/e^ll for an a xi symmetric disturbance 
in place of throe for the no n- a xi symmetric diS' 
turbance . 

We will, therefore, consider the two cases separately for 
satisfying the boundary conditions at the inner wall. 


Axlsymmetric Disturbances (n = 0) 

As explained above, the eigenfunctions for an 

axisymmetric disturbance can be expressed in terms of only 

two constants U,| and 5 for example, the axial velocity 

eigenfunction v (r) may be written as 

z 

v^(r) = ^1 ^1’ (3«26) 


The boundary conditions at the inner wall require 
that the disturbance velocity components be zero at r = Y , 
so tloat in view of equation (3.26), we have 


Vr1 ' V ) 

Vr2 ( Y ) 


'> 




0 : (3.27) 


If det 'stands for the determinant of the coefficient 
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matrix, it is required, that det he zero for a non-trivial 
Solution of the set of homogeneous oqua.tions (3*27). Since 
the determinant det is, in general, complex, we need to 
ensure that the absolute va,lue of det be zero, that is, 

dot I = 0. (3*23) 

Further, since the elements of the determinant in equation 
( 3 . 28 ) are functions of k,uj and R, only certain combina- 
tions of these parameters will allow this equation to be 
satisfied, thus yielding the eigenvalue k (for spatial 
analysis) or the eigenvalue U) (for temporal analysis). 

Once an eigenvalue is found, the eigenfunctions v^(r), 

V (r) and p(r) are easily calculated in terms of only one 
z 

arbitra.ry constant, say , since P.j and U..| are now rela- 
ted by equation (3.27) as 




( 3 . 29 ) 


For convenience, we set = (1 .0 + i 0.0), and find all 
the eigenfunctions from equations of the form given in 
equation (3*26). 


3.3.2 Fon-axisymmetric Disturbances (n 0) 

Following an analysis similar to that above for 
axisymmetric disturbances, it is obvious that the bound.ary 
conditions at the inner wall will bo satisfied if 
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Vj,2( i ) 

^r3 

VggY) 

VggCY ) 

■^03 


) 

^z3 


* 

= 0. (3.30) 


The determinant, x/liose absolute value should he zero, is, 
therefore, of order 3 in contrsist to that of order 2 in the 
case of an axisyiiuiietric disturbance. Hius, the complex 
eigenvalue is again found by satisfying equation (3.28) 
where det is now given by the following 





"^rl 

( y ) 

y ) 

b3^ V ) 



det 

= 

■^01 

( V) 

Yq2( Y) 

Yg3(y) 

(3.31) 




^z1 

( Y) 



• 

The 

ratios 

Of 

and 

. to 

are given by (see 

equation 

(3.30)) 










'i ) 

v^i(y) 


y^^CY) 





t) 


- ’^93'- ^ 


Y O *0^ J 

and 











r) 

b2^ If ) 

- YggC Y ) 






Y) 


- ^03^ ^ ^ 


V o i 


so that arbitrarily setting = (1 ►O + i 0.0),, the eigen- 
functions can be found from, equations of the form given 
in equation (3.21). 
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EESULTS AIID DISOJSSIOIT 

The disturbance equations were solved on an 
lEN 70^4 Computer using double precision mode. The com'- 
puter programrae, included in Appendix III 5 can be used 
for studying both the spatial and temiooral stability of 
flow in a concentric annulus to infinitesimal axis ymme trie 
and non-axisynmietric disturbances . However, results are 
presented here only for the spatial stability of annular 
flow to axisyrmiietric disturbances. The dimensionless 
frequency co iNias given three different values 0.1, 0.7 
and 1 .0 while the modified Reynolds number R vra.s varied 
from 1,000 to 10,000. Three different anmili were con- 
sidered with diameter ratio Y =0.2, 0.7 and 0.8. 

4,1 Comuutational Procedure 

It is an unfortunate fact that the values of 
the constants (V, U and P’s in equations (3.16), (3.18) 
and ( 3 . 19 )) increase as m,ore terms are retained in the 
series expansions of the eigenfunctions. Fortunately, 
however, the exponent of (1-r) also increases rapidly 
and since a term in the aeries consists of the product 
lihe ( l-r)^^^”*”' \ it was found that for (1-r) less 
than 1, the successive .terms in the series decrease in 
magnitude. Then if r^ represents the preassigned value 
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cf* r up .to which the series solution was carried start- 
ing from the outer wall, the criteria used for termina- 
ting the series was that for r = r , the rsitio of the 
last term retained to the partia.l sum up to a tern preced-^ ■ 
ing the last one must he less tiian a preassigned epsilon. 
This ratio was found for all the three series (equations 
(3*16), (3«1-8) and ( 3 . 19 )) and the maximum of these ratios 
was required to he less than epsilon. In view of the 
accuracy requirements, a value of 10 was used for 
epsilon in all cases. It was found that the nuniher of 
terms required in the series varied from 9 to 25 tor a 
tyiDical r„ of 0.9; the higher nuiiiher of terms required 
for a larger value of the parameter 11. Hie values of the 
basic flow velocity and its derivative as calculated h;/ 
the series (equation (3*1)) and hy the actual formula 
(equation (2.l4)) were compared and found to he the sane 
up to at least seven significant figures even if only 
five terms were retained in series. The step size used 
for the fourth order Runge-Kutta method v/as taken as 


0.005. 


It is a well known fact [8, 9 Jthat the error 
involved in the Runge-Kutta method is not easily trace- 
able. Collatz [s, page 71 ] has provided a rough guide 
for finding a reasonable length of the .inte^gratioii step 
but in application to this problem, it l^Jas found that 



32 


Ills criteria was too stiff to te practicable. Ralstonf9] 
has provided some error estimates for the Eunge-Kutta 
method but they are good only for a single differential 
equation while xre have three coupled differential equa- 
tions for an axis ymrae trie disturbance. It may be noted, 
however, that up to E = 10,000, no numerical instability 
was encountered in any of the cases studied. 


4-. 2 Stability 


Since the floir is spatiall^r unstable if the 
real part of the complex wave nmaber k is greater than 
zero, regions in the first quadrant of the complex k-plane 
were explored for any possible eigenvalues by the search 
technique described in Appendix II. Some of the regions 
examined are shown in figure V.1. No eigenvalue was 
found to exist inside these regions concluding, thereby, 
that for an infinitesimal disturbance propagating down- ‘ 
stream all the eigenvalues for R ^ 10,000 had negative 
real parts. Ihus, the flow in a concentric annulus was 
found to be spatially stable to infinitesimal axisyamie- 
tric disturbances up to a modified Reynolds nuiiiber of 
10,000. These results arc In agreement with those of 
Mott and Joseph for temporal stability 
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^♦3 The Least S talkie Mode 

It may be recalled that for a linear stability 
analysis, every component of the double Fourier series 
forming the arbitrary disturbance (cf. equation (1.7)) 
can be considered independently, Hius , even though 
there are an infinite number of modes governing the pro- 
pagation of the disturbance, the most predominant one is 
the least stable mode (one for -which | | is minimum) 

or an unstable mode (if one exists). Results are, there- 
fore, presented here (in figures 4*2 throu^ 4,10) for 
the least stable mode since no unstable mode -was found 
for ann-ular flow up to B = 10,000. 

Figures 4*2 through 4.4 show, the variation of 
k„ with B (or B^) for the least stable mode with u) as a 
parameter. Each figure corresponds to a different dia- 
meter ratio y , For a particular Y , it is seen that 
each - B curve has two branches . Below a certain 
Reynolds number for the least stable mode follows 

OIK curve and above this Reynolds number, it follows ano- 
ther curve. The Reynolds number, at which such a change 
of branch occurs, is itself a function of . As ui is 
increased, this Reynolds number shifts to a lower value. 
This behaviour is found for all the three diameter ratios 
except for a slight abnorimlity for y = 0.5. It is 
also to be noted that for a fixed frequency and diameter 
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ratio, the least stable mode becomes less stable (since 
j I decreases) as the Reynolds nmiber Increases. This 
suggests a possibility of there being a critical Reynolds 
number above uhich the flow may be unstable. Figures 
4,2 through 4.4 indicate, however, that even if it exists, 
the critical Reynolds number wouCLd be much larger than 
10,000 - the maximum TOlue of R considered here. These 
figures also shov/ that the anniilar flow is more stable 
to a higher frequency disturbance for a fixed Reynolds 
numljer and diameter ratio.] the effect (of u; on k^) beco- 
ming less predominant as the diaideter ratio increases. 
Also, the effect of increasing the diameter ratio of the 
a.nnulus is to make the flow more stable for a giveii u) and 
R, hott and .Joseph [6 ]also found a similar behaviour in 
the case of temporal stability. 

The variation of k^ with R (or R. ) for the 

1 d 

least stable mode is shown in figures 4.5 through 4.7 

f:fgw’e 

with u> as a parameter; eaclVcorrespondi ng to a different 

t 

diameter ratio V . These figures are quite similar to 
those discussed above. Hero also, for ea,ch ui and Y , 
the - R curve is found to have two branches. However, 
for a fixed frequeneji" and diame.tor ratio the mvolength 
(= 1/kj_) remains; almost constant as R increases ten fold 
while I k^ I decreases by a factor of about four for the 
same rise in R. Tlie dimensionless phase velocity 
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of the least stable mode is nearly unity 
(see tables 1 through ^-3) . It may be observed that 
for a particular oo and Y , increases with E and tends 
to a nearly constant value, Hiis increase is maximum 
around the Reynolds number t/here a break in the k^ - R 
curves occurs-. Moreover, for a fixed Reynolds nuraber, the 
phase velocity increases -^d-th W ( V held constant) and 
also with 'Y C held cons tant ) . 

Figures ^+-*8 through ^,10 show the variation'of 
k^ R with >U)R for different diameter ratios. These cur** 
ves were plotted to see whether; as in the case of pipe 
Poiseuille flow [^] , any simple relation exists between 
1{R and u>R-, If such a relation does exist, coR would be 
the governing parameter rather than the frequency oo and 
the modified Reynolds number R independently. Figures 
through ^.10, however, indicate that, contrary to 
that for flow in a pipe, u)R cannot be taken as the go- 
verning parameter for flow in a concentric annulus , 

These figures show that as the diameter ratio increases 
the number of branches of the k^R - tOR curve decreases. 
Some results are also tabulated in tables ^+-1 
through 4-3 , These tables display a few modes that are 
more stable than the corresponding least stable mode. 
Though some of these more stable modes are very dose to 



the corresponding least stable mode (cf. S .Kos . 9 and 12 
table 4-1 5 S.l'Ios. 3 and 16 table 4-25 nnd S ,1'Ios . 9 and 
l4 table 4 - 3)5 eigenvalne searcli teclinioue coupled 

uitii the Iterative scheme was able to isols tc them with- 
out much difficult 3 r. 
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Eigenvalues for the Least Stable Mode and for some Other 

Stable Modes ( Y =0.2) 
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CORCLUSIOl'lS 

The effect of modified Reynolds number R, 
diamotor ratio Y , and frequenc^r u> on the relative 
stability of concentric annular flow and on the wave- 
length of tlio disturbance can be summarized as follows ; 

(i) For a fixed Y and uj , the least stable mode 
becomes less stable and has a larger \\rave- 
length as R is increased. 

(ii) For a fixed Y and R, the least stable mode 
becomes less sta.ble and IjOS a larger*' wave- 
length as is decreased. 

(iii) For a fixed to and R, the least stable mode 
becomes loss stable a.nd j..as a larger wave- 
Icngtli as Y is decreased. 

On the basis of these results, it is concluided 
that up to a joodifiod Ilejmolds numbor of 10,000, the 
concentric annular flow is spatially stable to infini'ce- 
simal axis ynmic trie disturbances. These results are in 
conformity with those of hott and Joseph for the tempo- 
ral stability [6 j . 
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STEP-BY-STEP INTEGEATION TECHNIQUE 

Tlie fourth order Runge-Kutta method used for 
the step-hy-step integration of the coupled differential 
eauations (3.22) through ( 3 * 25 ) is given in a tabular form 
here. Tables AI-1 and AI-2 exhibit the appropriate for- 
mulae for first and second order differential equations. 

In these tables h is the stop size, D is ■ the operator 
d/dx, and x^, y^ and Dy^ are the values known at the start- 
ing point of the integration step. The procedure is to 
calculate first and then go to the first column in the 
second row. After calculating K 2 ty traversing the second 
row from left to right, we go to the first column in the 
third rov;, and so on till is calculated. The factors 
in the 'Correction' column are then calculated before we 
go to i'if th rov;; X-j , y^ and Dy^ are the values at the end 
of the present step, that is they are the starting point 
values for the next integration step. The whole procedure 
is, therefore, repeated for successive steps. 


l.f.T. KANPUR 

gENTR-AL I raRARY 
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Eunge-Kutta Scheme for Differential Equations of the 
First Order; Dy = f(x,y) 
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APPENDIX II 


EIGEWALUE SEARCH TECHNIQUE 

The eigenvalue search technique developed "by 
Gang and Rouleau [^4- J uses the result of CauGh 3 ?' Integial 
Theorem for determining the numher of zeros of a complex 
function in any closed region of the argument in which 
txiG function is analj^'tic* 

For a function f(z), analytic except for poles 
in the interior of a closed curve C on the z-plane, and 
for fCz) and I)f(z) continuous on C, the residue theorem 
gives 

dz , (AII.1) 

'.fliOTQ Z and P denote respectively the numher of zeros and 
tlie numher of poles of f(z) within the closed region C 
(counted vdtJi t.’icir multiplicities), and where D is the 
operator d/dz. Since Df(z)/f(z) is the derivative of 
In [fCz)] , to calculate the definite integral the 
right-hand side of equation (AII.1), it is sufficient to 
knoi/ only tlie variation of 

In I f(z) j + i angle [f(z) ] 

when tho variahle z describes the contour C in the 
positive sense, that is, in the counter-clockwise sense. 
But |f(z) I returns .to its initial value, while the angle 
of f(z) increases by 2 tt M, M being either zero or a 




positive or negative integer. Tims equation (AII.1) 
reduces to 


2 IT Ml _ 
2iTi 


(All. 2) 


that is, the difference (Z-P) is equal to the integer 
quotient ohtoined hj?- dividing the variation of the angle 
of f(z) by 2¥ as the variable z describes tie bounds igy 
C in the positive sense. 

The function f(z) can be expressed in terms of 
its real and iau ginar;/' parts as follov/s 

f(z) = X + iY. 

Then, as the point z = x + iy describes the curve C in 
the positive sense, tlie point vrlioso coordiiiates are 
(X, y) with respect to a S 3 .rsten of rectangiilar axes uith 
the same orientation as the first sj^sbem, describes also 
a closed curve C^, and it is only required to count the 
numl-or of revolutions v/hich the radius vector joining tlie 
origin of coordinates to the point (X, Y) has turned th- 
rough in one sense or the other. There is thus no need 
of plotting the actual curve . Instead, the integer M 
is given b 3 ^ the net multiples of 2 it by which the phase 
angle of f(z) changes in going round the curve . ten 
the point z describes the curve C in the positive 
sensCy anj'' counterclockv/isc movement of the phasor for 
f(z) is considered positive while the cloclcwise movement 
Is considered negative. 



Zn s.ppliCcL'fcion to this prohiciiij Gsich 6 ioi"iG?"'.t 
of tlie determinant (equation (3.31)) is a function of 
the complex eigenvalue k or u . In fact, the determinant 
det can he considered as some high order poljrnoniial in 
the eigenvalue due to the recurrence relations for the 
coefficients i:-:i the series solution of the stability equa- 
tions . Therefore, the function det (k) or det (w) has 
no poles, tliat is, P = 0 in equations (AII0I) and (All . 2 ), 
and Z = K. Thus, the prohlem of finding the number of 
eigenvalues ijithin a closed region of the complex k or to- 
piano is as sMplo as counting the nnt multiples of 2 Tr 
by uhich the phase angle of the determinant changes as k 
or ui assumes values on tJie closed contour in the complex 
plane . 

This method can also be used to provide a close 
apnroxiriation to the true eigenvalue so tint an iterative 
tcchniauc may converge to it. If it so happens that the 
closed contour chosen on the complex k or to -plane passes 
close to an eigenvalue within the contour at some location, 
the phase angle of the determinant \.dll change hj about 
180^ and the magnitude of the determinant idll show a 
dip. Ilovrever, depending on the closeness of the boundarj;^ 
to the eigenvalue, this litiy or may not be the case, never- 
theless, ve can almys subdivide the region investigated 
on the complex plane and thus isolate a relatively small 
region which contains tiie eigenvalue . Experience has 
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shown that except for the least stable mode, it is 
sufficient to know only significant digits correctly 
in the approxanate eigenvalue. The secant method used 
for the iterative minimization scheme then converges to 
the true eigem'alue easily. For the least stable mode, 
it is generally required to know the approximate eigen- 
TOlue to one noro order of accuracy. 



double precision form. Although tie double precision 
complex arithmetic is not available on the IBII 7044 
s^'-stera it uas created out of necessity. Thus, every 
double precision complex variable was expressed in terns 
of tv;o double i-srecision variables w'^ich represent the 
real .and i.. a^;iuar 3 r parts of the complex que..ntity. 



G- ^ Uf\j 



